Strange Quark Contributions to Parity- Violating Asymmetries in the Forward GO 

Electron-Proton Scattering Experiment 



in 



D. S. Armstrong, 1 J. Arvieux, 2 R. Asaturyan, 3 T. Averett, 1 S. L. Bailey, 1 G. Batigne, 4 D. H. Beck, 5 
E. J. Beise, 6 J. Benesch, 7 L. Bimbot, 2 J. Birchall, 8 A. Biselli, 9 P. Bosted, 7 E. Boukobza, 2 ' 7 H. Breuer, 6 
R. Carlini, 7 R. Carr, 10 N. Chant, 6 Y.-C. Chao, 7 S. Chattopadhyay, 7 R. Clark, 9 S. Covrig, 10 A. Cowley, 6 
D. Dale, 11 C. Davis, 12 W. Falk, 8 J. M. Finn, 1 T. Forest, 13 G. Franklin, 9 C. Furget, 4 D. Gaskell, 7 J. Grames, 7 

K. A. Griffioen, 1 K. Grimm, 1 - 4 B. Guillon, 4 H. Guler, 2 L. Hannelius, 10 R. Hasty, 5 A. Hawthorne Allen, 14 
T. Horn, 6 K. Johnston, 13 M. Jones, 7 P. Kammel, 5 R. Kazimi, 7 P. M. King, 6 ' 5 A. Kolarkar, 11 E. Korkmaz, 15 
W. Korsch, 11 S. Kox, 4 J. Kuhn, 9 J. Lachniet, 9 L. Lee, 8 J. Lenoble, 2 E. Liatard, 4 J. Liu, 6 B. Loupias, 2,7 
A. Lung, 7 G. A. MacLachlan, 16 D. Marchand, 2 J. W. Martin, 10 - 17 K. W. McFarlane, 18 D. W. McKee, 16 
R. D. McKeown, 10 F. Merchez, 4 H. Mkrtchyan, 3 B. Moffit, 1 M. Morlet, 2 I. Nakagawa, 11 K. Nakahara, 5 M. Nakos, 16 

R. Neveling, 5 S. Niccolai, 2 S. Ong, 2 S. Page, 8 V. Papavassiliou, 16 S. F. Pate, 16 S. K. Phillips, 1 M. L. Pitt, 14 
M. Poelker, 7 T. A. Porcelli, 15 ' 8 G. Quemener, 4 B. Quinn, 9 W. D. Ramsay, 8 A. W. Rauf, 8 J.-S. Real, 4 J. Roche, 7 ' 1 
£sj ! P. Roos, 6 G. A. Rutledge, 8 J. Secrest, 1 N. Simicevic, 13 G. R. Smith, 7 D. T. Spayde, 5 - 19 S. Stepanyan, 3 

M. Stutzman, 7 V. Sulkosky, 1 V. Tadevosyan, 3 R. Tieulent, 4 J. van de Wiele, 2 W. van Oers, 8 E. Voutier, 4 
W. Vulcan, 7 G. Warren, 7 S. P. Wells, 13 S. E. Williamson, 5 S. A. Wood, 7 C. Yan, 7 J. Yun, 14 and V. Zeps 11 

(GO Collaboration) 

'Department of Physics, College of William and Mary, Williamsburg, VA 23187 USA 
2 Institut de Physique Nucleaire d'Orsay, Universite Paris-Sud, F-91406 Orsay Cedex FRANCE 
• 3 Yerevan Physics Institute, Yerevan 375036 ARMENIA 

■* Laboratoire de Physique Subatomique et de Cosmologie, F-38026 Grenoble Cedex FRANCE 
I , 5 Loomis Laboratory of Physics, University of Illinois, Urbana, IL 61801 USA 

■ Physics Department, University of Maryland, College Park, MD 20472 USA 

' 'Thomas Jefferson National Accelerator Facility, Newport News, VA 23606 USA 

\ 8 Department of Physics, University of Manitoba, Winnipeg, MB R3T 2N2 CANADA 

■ 9 Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213 USA 

' 10 Kellogg Radiation Laboratory, California Institute of Technology, Pasadena, CA 91125 USA 

11 Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506 USA 
'■ 12 TRIUMF, Vancouver, BC V6T 2A3 CANADA 

jjj ■ 13 Department of Physics, Louisiana Tech University, Ruston, LA 71272 USA 

I ' 1J> Department of Physics, Virginia Tech, Blacksburg, VA 24061 USA 

15 Department of Physics, University of Northern British Columbia, Prince George, BC V2N 4Z9 CANADA 
_. 1 16 Physics Department, New Mexico State University, Las Cruces, NM 88003 USA 

11 Department of Physics, University of Winnipeg, Winnipeg, MB R3B 2E9 CANADA 
18 Department of Physics, Hampton University, Hampton, VA 23668 USA 
19 Department of Physics, Grinnell College, Gnnnell, IA 50112 USA 
(Dated: February 3, 2008) 

We have measured parity-violating asymmetries in elastic electron-proton scattering over the 
range of momentum transfers 0.12 < Q 2 < 1.0 GeV 2 . These asymmetries, arising from interference 
of the electromagnetic and neutral weak interactions, are sensitive to strange quark contributions 
to the currents of the proton. The measurements were made at JLab using a toroidal spectrom- 
eter to detect the recoiling protons from a liquid hydrogen target. The results indicate non-zero, 
Q 2 dependent, strange quark contributions and provide new information beyond that obtained in 
previous experiments. 

PACS numbers: 11.30.Re, 13.60.-r, 14.20.Dh, 25.30.Bf 



At short distance scales, bound systems of quarks have 
relatively simple properties and QCD is successfully de- 
scribed by perturbation theory. On size scales similar 
to that of the bound state itself, ~ 1 fm, however, the 
QCD coupling constant is large and the effects of the 
color fields cannot yet be calculated accurately, even in 
lattice QCD. In addition to valence quarks, e.g., uud for 
the proton, there is a sea of gluons and qq pairs that 
plays an important role at these distance scales. From 



a series of experiments measuring the neutral weak scat- 
tering of electrons from protons and neutrons, we can ex- 
tract the contributions of strange quarks to the ground 
state charge and magnetization distributions (e.g., mag- 
netic moment) of the nucleon. These strange quark con- 
tributions must originate in fluctuations of gluons to ss 
pairs because there are no strange valence quarks in the 
nucleon. There have been numerous estimates of strange 
quark contributions to nucleon properties within various 
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phenomenological models and also in state-of-the-art lat- 
tice calculations 0, ; many focus on the contribution 
to the magnetic moment. In this paper, we report on 
a new measurement sensitive to strange quark contribu- 
tions over a range of distance scales. 

Separation of the strange quark contributions to nu- 
cleon currents in the context of the neutral weak interac- 
tion dates back to Cahn and Gilman Q and was devel- 
oped by Kaplan and Manohar Q. Because the coupling 
of both photons and Z bosons to point-like quarks is 
well defined, it is possible, by comparing the correspond- 
ing currents^to separate the contributions of the various 
flavors [ilia,!]]]- The charge and magnetic form factors of 
the proton can be written [i = 7, Z) 

G P E,M — e% ' U G % E,M + e% ' {pE,M + ^E,m) > U) 

neglecting the very small contribution from heavier fla- 
vors. For the ordinary electromagnetic form factors the 
charges are e 7 = +2/3, —1/3 for u and d/s quarks, re- 
spectively. Assuming that the proton and neutron are re- 



where 

r=^, £ =(l + 2(l + r)tan^) \ 

V = e{G E f + r(G] f ) 2 , and s' - y/r(l + r)(l - e 2 ), 

Q 2 is the squared four- momentum transfer (Q 2 > 0), Gf 
and a the usual weak and electromagnetic couplings, M p 
the proton mass and 9 the laboratory electron scattering 
angle. The three new form factors in this asymmetry, 
G§, Gfj and G e A may be separated by measuring elas- 
tic scattering from the proton at forward and backward 
angles, and quasi-elastic scattering from the deuteron at 
backward angles 0. 

The GO experiment was performed in Hall C at 
Jefferson Lab. We used a 40 fiA polarized electron 
beam with an energy of 3.031 ± 0.001 GeV over the 
measurement period of 700 h. It was generated with 
a strained GaAs polarized source 0] with 32 ns pulse 
timing (rather than the standard 2 ns) to allow for time- 
of-flight (t.o.f.) measurements. The average beam po- 
larization, measured with a M0ller polarimeter in 
interleaved runs, was 73.7 ± 1.0%. Helicity-correlated 
current and position changes were corrected with active 
feedback to levels of about 0.3 parts-per-million (ppm) 
and 10 nm, respectively. Corrections to the measured 
asymmetry were applied via linear regression for residual 



lated by a simple exchange of u and d quarks |8| (and the 
corresponding anti-quarks), the ordinary neutron form 
factors can be written in terms of these same contribu- 
tions 

2 1 

r^ n n — r< d (r m _i_ r< s \ fo\ 

U E,M — ^ U E,M _ g \ U E,M ' U E,MJ • \ z ) 

A complete separation of the G q E M , and, in particular, 
isolation of G S E M , requires a third combination. In this 
paper, new measurements of the weak interaction form 
factors of the proton are presented which allow us to 
determine the strange quark contributions. These form 
factors are written (Eqn.^) in terms of the weak charges, 
e z = 1 — 8/3 sin 2 0vK 7 —1 + 4/3 sin 2 6>w for the u and d/s 
quarks, respectively, where 9\y is the weak mixing angle. 

In order to isolate the small contribution to elastic 
electron-proton scattering from the neutral weak current, 
we measure the parity-violating asymmetry for longitu- 
dinally polarized (i? and L) electrons 
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helicity-correlated beam current, position, angle and en- 
ergy variations and amounted to a negligible total of 0.02 
ppm; the largest correction was 0.01 ppm for helicity- 
correlated current variation. We made one further cor- 
rection of, on average, +0.71 ± 0.14 ppm to the asym- 
metries in all detectors (~ 5% variation from detector to 
detector). It was associated with a small (~ 10 -3 ) frac- 
tion of the beam current with a 2 ns structure ( "leakage 
beam" : tails of beams from other operating halls) and a 
large charge asymmetry (~ 570 ppm); it was measured 
in otherwise 'forbidden' regions of the t.o.f. spectra. 

The polarized electrons scattered from a 20 cm liquid 
hydrogen target 01; the recoiling elastic protons were 
detected to allow simultaneous measurement of the wide 
range of momentum transfer, 0.12 < Q 2 < 1.0 GeV 2 . 
This was effected using a novel toroidal spectrometer de- 
signed to measure the entire range with a single field 
setting and with precision comparable to previous exper- 
iments. The spectrometer included an eight-coil super- 
conducting magnet and eight sets of scintillator detec- 
tors. Each set consisted of 16 scintillator pairs used in 
coincidence to cover the range of momentum transfers 
(smallest detector number corresponding to the lowest 
momentum transfer). Because of the correlation between 
the momentum and scattering angle of the elastic protons 
(higher momentum corresponds to more forward proton 
scattering angles), detector 15 covered the range of mo- 
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FIG. 1: Example of the raw asymmetry, A mea s, (data points) 
and yield (histogram) as a function of t.o.f. for detector 8. 

mentum transfers between 0.44 and 0.88 GeV 2 which we 
divided into three t.o.f. bins with average momentum 
transfers of 0.51, 0.63 and 0.79 GeV 2 . For the same rea- 
son detector 14 had two elastic peaks separated in t.o.f. 
with momentum transfers of 0.41 and 1.0 GeV 2 ; detector 
16, used to determine backgrounds, had no elastic accep- 
tance. Custom time-encoding electronics sorted detector 
events by t.o.f.; elastic protons arrived about 20 ns after 
the passage of the electron bunch through the target (see 
Fig. The spectrometer held integral and ultimately 
the Q 2 calibration (AQ 2 /Q 2 = 1%) was hne-tuned using 
the measured t.o.f. difference between pions and elastic 
protons for each detector. All rates were corrected for 
dead-times of 10—15% on the basis of the measured yield 
dependence on beam current; the corresponding uncer- 
tainty in the asymmetry is ~ 0.05 ppm. Standard radia- 
tive corrections 01 nl the range of 1 - 3%, determined 
by comparing simulations with and without radiation, 
were also applied to the asymmetries. Lastly, there is 
an uncertainty of 0.01 ppm due to a small component of 
transverse polarization in the beam. 

As shown in Fig. ^ a background extends on both 
sides of the elastic proton peak at a t.o.f. ~ 20 ns. This 
background is essentially all protons (as determined from 
energy loss measurements in a sampled data set): quasi- 
elastic protons from the aluminum target windows and 
inelastic protons from both the hydrogen and the alu- 
minum. The measured asymmetry has two components 

A m eas = (1 ~ f) A el + f A back (4) 

where A e i is the raw elastic asymmetry and / is the back- 
ground fraction; in the actual analysis t.o.f. fits to the 
yield and asymmetry in the region of the elastic peak 
are used. The yield is typically modeled with a Gaussian 
elastic peak and a polynomial background. The asym- 
metry model comprises a quadratic background and a 
constant for the elastic. For higher numbered detectors 



TABLE I: Asymmetries and uncertainties measured in the 
present experiment p^l . The contributions to the systematic 
uncertainties are summarized in Table HTl 
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TABLE II: Systematic uncertainties for measured asymme- 
tries. The first six uncertainties are global, deadtime is point- 
to-point and the background is a combination (see text). 



Source Uncertainty 

Helicity-correlated beam parameters 0.01 ppm 

Leakage beam 0.14 ppm 
Beam polarization 1.0% 
Ordinary radiative corrections 0.3% 

Transverse polarization 0.01 ppm 
Q 2 1% 

Background correction 0.2 - 9 ppm 

Deadtime 0.05 ppm 



the background asymmetry is positive. In particular, for 
detector 15 the background asymmetry has a maximum 
value of about 45 ppm in the region of the elastic peak. 
As substantiated by a Monte Carlo simulation, this pos- 
itive asymmetry is caused by a small number of A and 
£ weak-decay protons scattered inside the spectrometer 
magnet. The smooth variation of the region of positive 
asymmetries is tracked from detectors 12-14 through to 
detector 16; the background asymmetry for the large ac- 
ceptance of detector 15 is then corrected by interpolating 
these measured background asymmetries. As a check, 
the same fitting procedure described above is also used 
for detector 15 and gives consistent results. 

The elastic asymmetries for the experiment, A p h ys {A e i 
corrected for all effects described earlier) are presented in 
Table U The statistical uncertainties include those from 
the measured and the background asymmetries. The sys- 
tematic uncertainties (Table ITTf) are dominated by those 
from the background correction. This uncertainty is es- 
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timated from the range of elastic asymmetries generated 
from a variety of different background yield and asymme- 
try models. These models are bounded by the measured 
slopes of background yields and asymmetries on cither 
side of the elastic peak and varied continuously between 
these limits. The uncertainty in the background asym- 
metry for detector 15 is conservatively taken to be the 
difference between interpolated background asymmetries 
in successive detectors as described above. We have also 
estimated the global and point-to-point contributions to 
these uncertainties from the extent to which a change in, 
e.g., the background asymmetry functional form, consis- 
tently changes the asymmetries in all the affected detec- 
tors. 

The results of the experiment are shown as a function 
of momentum transfer in Fig. [2J The quantity 
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(A p hys — Anvs) > (5) 



(where r\ (Q 2 ) = tG p m / sG p E ) is determined from the dif- 
ference between the experimental asymmetry and the 
"no-vector-strange" asymmetry, A^vs- Anvs is calcu- 
lated from Eqn.|3]with G% = G S M — for all values of Q\ 
and using the electromagnetic form factors of Kelly [15J . 
Also shown is the excellent agreement with the HAPPEX 
measurements Il7j made at nearly the same kinematic 
points (with small corrections to the asymmetries, < 0.2 
ppm, to adjust them to the GO beam energy). The error 
bars include the statistical uncertainty (inner) and statis- 
tical plus point-to-point systematic uncertainties added 
in quadrature (outer). The error bands represent, for the 
GO experiment, the global systematic uncertainties: from 
the measurement (upper) and from the uncertainties in 
the quantities entering A^vs (lower). These quantities 
are: the calculated value of the axial-vector form factor 
normalization |l8j (differing from gA/gv by electroweak 
radiative corrections), the same dipole momentum trans- 
fer dependence for G A (Q 2 ) as is deduced for Ga(GP) [13- 
the axial vector strangeness contribution As [20j, and 
the electroweak radiative corrections |2l|. The sensitiv- 
ity of the result to electromagnetic form factors is shown 
separately by the lines on the plot. For the alternative 
form factor parameterizations of Friedrich and Walcher 
(FW) [12 (dashed) and the combination (dotted): Ar- 
rington "Rosenbluth" [23| - proton, and Kelly 0] - neu- 
tron, the effective A^vs is shown (e.g., for the FW pa- 
rameterization, the value of G S E + r)G s M at Q 2 = 0.63 
GeV 2 increases from 0.059 to 0.072). Alternately, the 
uncertainties in the Kelly form factor fits would increase 
the width of the uncertainty band (lower) for Amvs a t 
each Q 2 by about 25% if included there. 

The G S E + r]G s M data shown in Fig. have a system- 
atic and intriguing Q 2 dependence. For reference we note 
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FIG. 2: The combination G S E + rjG s M for the present mea- 
surement. The gray bands indicate systematic uncertainties 
(to be added in quadrature); the lines correspond to different 
electromagnetic nucleon form factor models (see text). 



terize our result with a single number, we tested the hy- 
pothesis G S E + r)G s M — by generating randomized data 
sets with this constraint, distributed according to our 
statistical and systematic uncertainties (including corre- 
lated uncertainties). The fraction of these with \ 2 larger 
than that of our data set was 11%, so we conclude that 
the non-strange hypothesis is disfavored with 89% con- 
fidence. More important is the Q 2 dependence of the 
data. The initial rise from zero to about 0.05 is consis- 
tent with the findingthat G S M {CP^ = 0.1 GeV 2 ) ~ +0.5 
from the SAMPLE 24], PVA4 Q and HAPPEX [HI 
measurements. Because 7/ increases linearly throughout, 
the apparent decline of the data in the intermediate re- 
gion up to Q 2 ~ 0.3 indicates that G S E may be negative 
in this range. There is also some support for this conclu- 
sion from the combination of GO and PVA4 j2|| results 
at Q 2 — 0.23 GeV 2 . There is a significant trend, consis- 
tent with HAPPEX to positive values of G" E + 7]G S M 
at higher Q 2 . Experiments planned for Jefferson Lab, 
including GO measurements at backward angles, and 
MAM1 (Mainz) will provide precise separations of G S E 
and G S M over a range of Q 2 to address this situation. 

In summary, we have measured forward angle parity- 
violating asymmetries in elastic electron-proton scatter- 
ing over a range of momentum transfers from 0.12 to 1.0 
GeV 2 . These asymmetries determine the neutral weak 
interaction analogs of the ordinary charge and magneti- 
zation form factors of the proton. From the asymmetries 
we have determined combinations of the strange quark 
contributions to these form factors, G S E + r]G s M , which, 
together with other experiments, indicate that both G S M 
and G S E are non-zero. 
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